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This paper shows that the energetics of Boussinesq and anelastic fluids possesses a 
term that can be identified as the approximation SWb a to the compressible work of ex- 
pansion/contraction SW = —Pdv, where P is the pressure and v is the specific volume. 
It follows that Boussinesq and anelastic fluids admit explicit compressible effects and 
conversions between internal energy and mechanical energy, under the form of apparent 
changes in gravitational potential energy resulting from changes in density by diabatic 
and adiabatic effects. From the knowledge of 8Wb ai the corresponding approximation to 
the "heat" SQi, a can be constructed in a consistent way by requiring that the Maxwell 
relationships be satisfied, ultimately leading to the construction of a well defined approx- 
imation to the internal energy and ultimately of the full range of known thermodynamic 
potentials. These properties make it possible to endow common forms of the Boussinesq 
and anelastic approximations with fully consistent energetics and thermodynamics, even 
when diabatic effects and an arbitrary nonlinear equation of state for a binary fluid are 
retained, without loss of accuracy. In that case, it can be shown that the sum of kinetic 
energy and enthalpy is a conservative quantity, which plays the role of the total energy 
in the Boussinesq and anelastic approximations for both diabatic and adiabatic motions. 
This implies that gravitational potential energy can be regarded as the difference between 
enthalpy and internal energy, and hence as a pure thermodynamic property of the fluid. 

An importan t implication of the present results is to support the recent suggestion by 
iTailleuxl (|2009l ) that the Boussinesq approximation is capable of describing potentially 
large conversions between internal energy and gravitational potential energy in turbulent 
stratified fluids, which physically seems to require that an active role be played by the 
emission and dissipation of acoustic waves by molecular diffusive heating and cooling, 
as well as by the associated divergent velocity field, a surprising result with potentially 
important implications for our understanding of turbulent mixing in stratified fluids if 
further confirmed. Another implication of the results is to suggest that the form of the 
Boussinesq primitive equations currently used as the basis for a majority of numerical 
ocean models possesses a potentially significant spurious source of momentum, which can 
in principle be corrected by using an alternative and more physically-based definition of 
buoyancy in the hydrostatic approximation. 
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1. Introduction 



Many fluid flows of interest in engineering and geophysical fluid dynamics have low 
Mach number (M = Uq/c s ) and small relative density variations (p — po)/po, where Uq is 
a typical velocity scale, c s is the speed of sound, p is the density and pa a reference density. 
It has been common practice over the past century to regard such flows as incompressible 
or weakly compressible, where an incompressible fluid is one whose density dependence 
upon pressure is eliminated while still possibly retaining i ts dependenc e upon tempera- 
ture (and chemical composition as the case may be), e.g.. lLillv l (|l996h . There has been 
much effort in seeking to take advantage of the smallness of the two above parameters to 
develop sound-proof reduced sets of equations that are somehow simpler to study than 
the fully compressible Navier-Stokes equations. Two particular classes of approximations 
have been particularly influential and key to simplifying the numerical and theoretical 
analysis of low Mach number fluid flows, and will be under focus in thi s pap er. The first 
one is the Oberbeck-Boussinesq approximation (after lOberbeckl (Il879h and lBoussinesq 



(190i)), which in its most common form retains only the rotational divergence free com- 
ponent of the velocity field, and treats the density as constant everywhere except where 
it m ultiplies the accele r ation of gravity. The second one is the Anelast ic approximation , 



e.g.JOgura fc Phillips I (ll962h;lLipps fc Hemler I (ll982l):lBannon I (ll996l);|Purran I (Il989h: 
1982): Ir 



Ingersoll fc Pollard! (|l982F r ilngersoll I (|2005l ); iPauluis I (|2008l ). Many other sets can be 



constructed, which are beyond the scope of this paper, for insta nce by u s ing low Mach 
number asymptotics and mul t i-scal e expansion techniques, e.g., iMiiller I (Il998l ): iKlein I 
( 2009Ll2ni0h . lDavies fc etahl (|20f)3h offers a review of a number of commonly employed 
reduced sound-proof sets of equations, and show how they respectively represent normal 
modes on the sphere. 

The main focus of this paper is on how compressible effects and the coupling between 
mechanical energy and internal energy (the dynamics/thermodynamics coupling) are rep- 
resented in the Boussinesq/anelastic approximations, which are known to decouple either 
fully or almost fully the thermody namics from the dynamic s for adiabatic motions an d 
a linearised equations of state, i.e.. ISpiegel fc Veronis I (|l960l ): lOgura fc Phillips! (|l962h . 
In the latter case, the Boussinesq/ anelas t ic app roximations usually admit a well-defined 
conservative energy quantity, e.g.. iLillyl ( 19961 ). but the issue of the energetic and ther- 
modynamic consistency appears to be much less understood when diabatic effects and/or 
a nonlinear equation of state are retained, because the thermodynamics/dynamics cou- 
pling then becomes less trivial. In this respect, the oceanographic case is instructive, 
as oceanographers have used the Boussinesq hydrostatic approximation in conjunction 
with a realistic nonlinear equation of state (including the pressure dependence) for many 
decades as the basis for numerical ocean general circulation models of the kind used 
in climate studies, without any apparent obvious drawb acks apart from the lack of a 
well-defined and closed energy budget, see lTailleuxl(l2010h . Since Boussinesq ocean mod- 
els appear to work well with a "compressible" equation of state, one may wonder how 
essential the assumptions of "incompressibility" or "weak compressibility" are in the 
construction of the Boussinesq and anelastic approximations in the first place. 

To clarify this issue, a thermodynamics perspective is useful. Indeed, from a thermody- 
namic viewpoint, the very idea that it might be possible to simplify a particular equation 
of state to eliminate its dependence upon pressure while retaining its dependence upon 
temperature appears to become very dubious if not outright physically meaningless when 
diabatic effects are retained. This is because the way changes in density are affected by 
pressure changes depends critically on the particular thermodynamic transfor mation un- 
dergone by the fluid parcels. Take an isothermal transformation for instance. I Tailleux"! 



Energetics and thermodynamics of the Boussinesq/ Anelastic approximations 3 
(l2010l ) shows that the isothermal compressibility 7 = p 1 dp/dP\T is linked to the adia- 



batic compressibility p 1 dp/dP\ n = l/(pc 2 s ) (his Eq. A. 6) by: 
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where T = aT / (pc p ) is the adiabatic lapse rate. The key point here is that setting up the 
adiabatic compressibility to zero (by taking the zero Mach number limit c s = +00), while 
it makes the fluid effectively "incompressible" for adiabatic motions, fails in general to 
do so for nondiabatic transformations if the thermal expansion coefficient a is allowed 
to remain nonzero. Note that for seawater, typical values are: l/(pc 2 ) ~ 4.10 _10 P a _1 , 
whereas ar = a 2 T/{pc p ) « 7.5 x 10" 13 P a ~\ using c s = 1500 m.s -1 , a = 1Q- 4 K _1 , c p = 
4.10 3 J.K .kg -1 , and p = 10 3 kg.m~ 3 . For these values, the limit c s — +00 decreases 
the isothermal compressibility by about two to three orders of magnitude, so that even 
if it fails to fully eliminate compressibility effects for non-adiabatic motions, it appears 
nevertheless capable of reducing them considerably. This being said, if one agrees that the 
very concept of an "incompressible" fluid becomes physically meaningless when diabatic 
effects are retained, then one may also agree that it might not be that essential to 
set up the adiabatic compressibility to zero in the first place. Historically, the latter 
approach was originally motivated as a natural way to filter out sound waves, but it is 
now recognised that imposing c s = +00 is not necessarily to filter out sound waves, as 
such a filtering can be more simply achieved by usin g a hydrostatically adjusted pr essure 
in the equation of state for density, as discussed by Ide Szoeke fc Samelson I (|2002l ) . 

Whether Boussinesq/anelastic models accurately represent the conversions between 
mechanical energy and internal energy recently came to attention in relation with the 
question of how strong would the oceanic overturning circulation and meridional heat 
transport be, and how much diapycnal mixing would then be supported, if it were possible 
somehow to suppress the mechanical stirring due to t he wind and tide s, which ha s been a 
cont roversial issue f or the past decade, as reviewed in iTailleux I (jioOflh . iTailleuxl (|2010bl ) 
and iHughes fe al. I (120091) . The resulting configuration is often referred to as horizontal 
convection, see Hughes fc Griffiths (2008) for a recent review on the subject. The reason 
why understanding the nature of the thermodynamics/dynamics coupling is important 
in that case is because the steady-state mechanical energy balance reduces to: 

l P TSi Am = l^ < L2) 

which shows that in order to estimate the overall viscous dissipation rate, one has to 
estimate the overall work of expansion/contraction, where P is the pressure, v is the 
specific volume, dm = pdV is the elementary mass element of a fluid parcel, and pex is 
the viscous dissipation rate. As discussed by ITailleux I ( 2010b), the issue of estimatin 



B for a fully compressible fluid is a subtle one. As shown bv iPaparella fe Young I (|200 



and others, considerable analytical progress can be achieved for a Boussinesq fluid with 
a linear equation of state, as in that case, Eq. ()1.2|) becomes: 
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which can be integrated analytically, where k is the molecular diffusivity, go the accelera- 
tion of gravity and the terms within brackets are the surface area integra ted density at the 



bottom and top of the oceans respectively. Using typical oceanic values. I Wang fc Huang 



(2005) estimated the right-hand side to be 0(15 GW), which is at least two orders of mag- 
nitude smaller than the mechanical power input due to the wind and tides for instance 
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Figure 1. Experimental set-up suggesting that compressibility effects increase with degree 
of turbulence in stably stratified fluids. An Erlenmeyer flask is filled up with water at room 
temperature. A vertical stratification is then created by bringing the upper part of the water 
near boiling point, and a magnetic stirrer is dropped in the fluid. The top of the flask is then 
sealed up with a cork through which a thin tube is inserted in order to magnify the variations 
of the air- water interface. In laminar conditions (panel a), the fluid is near resting conditions, 
and the air-water interface moves downward very slowly as the result of molecular diffusion 
and/or radiative cooling. Activating the magnetic stirrer generates strong turbulence (panel b) 
that results in the spectacular drop of the air-water interface, presumably as the result of the 
contraction upon mixing discussed in the main text. This experiment was demonstrated to the 
author by Peter Rhines during a visit in Seattle in June 2008. 



The smallness of this number have been widely interpreted as evidence that the surface 
buoyancy fluxes could not by themselves be responsible for an overturning circulation 
and associated meridional heat transport of the observed magnitude. 

The smallness of the r.h.s. of Eq. (11.31) had been previously interpreted in the context 
of stratified turbulence (following I Winters fc al ( 1995t )) as implying that conversions be- 
tween internal energy and m echanical energy enter the energetics of turbulent stratified 
mixing only at second order. iTailleux I ([20091 ) argued, however, that the apparent small- 
ness of the thermodynamic/coupling in Boussinesq stratified turbulence actually hides 
two large and opposite conversions between internal energy (IE) and gravitational poten- 
tial energy (GPE) that almost exactly cancel out, viz., one conversion transferring IE into 
background GPE, associated with the smoothing out of the vertical reference tempera- 
ture gradient and a significant overall volume reduction (assuming that a increases with 
temperature, as is the case for water), the other conversion dissipating available GPE 
into IE in a way analogous to the viscous dissipation of kinetic energy (KE) into IE, and 
associated with the increase of the mean thermodynamic equilibrium temperature of the 
fluid but with only a negligible overall volume expansion. These results suggest therefore 
that compressibility effects must increase with the degree of stratification and turbu- 
lence, which appears to be supported empirically by the kind of laboratory experiment 
illustrat ed in Fig. [H but have yet t o be widely ac c epted , as most subseque nt studies so 
far, e.g., IWinters fe Young I (|2009T ): iHughes fe all (|2009t l ; iNvcander I (|2010h continue to 



favour the classical view that internal energy and thermodynamics/dynamics coupling 
play only a minor role in the energetics of turbulent stratified mixing. 
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At least two important misconceptions about weakly compressible fluids appear to 
greatly confuse the understanding of the role of compressible effects in stratified turbu- 
lent fluids. The first misconception is associated with the widespread (erroneous) belief 
that a weak thermodynamics/dynamics coupling is an intrinsic feature of all weakly com- 
pressible fluids irrespective of the particular thermodynamic transformations undergone 
by fluid parcels, whereas as far as we can judge, such a weak coupling is an intrinsic 
feature of adiabatic motions only, with no physical basis for this to be the case for di- 
abatic motions. The second misconception is associated with the widespread confusion 
about what a "weakly compressible" or "incompressible" fluid is actually supposed to 
be. Should such a property pertain to a fluid whose density dependence o n pressure is 
eliminated, as seems to have been the original intention, e.g.. iLillv I (|1996l ). or with the 
use of the constraints V • v = or V • (pov) — 0, as seems to be increasingly assumed 
following the realization that Boussinesq and anelastic fluids can be used with nonlinear 
equations of state? Based on the present analysis, the idea of an incompressible fluid 
seems to be justified only for a fluid with zero adiabatic compressibility in the context 
of purely adiabatic motions, but physically meaningless otherwise. Physically, the con- 
straints V-v = and V-(po v ) — do not imply incompressibility or weak compressibility, 
in contrast to what is often assumed, because they do not preclude diabatic or adiabatic 
density changes along fluid trajectories. The new view that emerges from the above argu- 
ments is that the Boussinesq and anelastic approximations are actually able to support 
potentially large compressible effects and conversions between internal energy and me- 
chanical energy. If so, this raises many questions that the present paper seeks to clarify. 
1) How do the conversions between internal energy and mechanical energy in a Boussi- 
nesq/ anelastic fluid compare with that of a fully compressible fluid? Does the answer 
depend on whether diabatic effects and/or a nonlinear equation of state are retained? 2) 
Do the Boussinesq and anelastic approximations conform with the first and second laws 
of thermodynamics? What is the form of the thermodynamic potentials supported by 
such approximations? How different are they from their exact counterparts? 3) Is it an 
intrinsic problem that many Boussinesq and anelastic models fail to be energetically and 
thermodynamically consistent, or can such models be modified to correct the problem? 
In the latter case, can this be done without modifying the formal order of accuracy of 
the original Boussinesq/ anelastic approximations? 



The main purpose of this paper is to provide an answer to all above questions, by 
showing that the Boussinesq and anelastic approximations can be endowed with well 
defined energetics and thermodynamics that closely mimic that of the fully compressible 
Navier-Stok es equat i ons. A couple of recent papers have touched up on some of thes e 
issues. Thus. I Young I (|2010l ) showed, using ideas previously developed by llngersoll I (|2005j ). 
that the seawater Boussinesq equation have a well-defined conserved energy quantity for 
an arbitrary no nlinear equation of state, although whether this extends to the diabatic 
case is not clear. Pauluis I (|2008t ) addressed a similar issue for the anelastic approximation 
for a binary fluid such as moist air, also with a highly nonlinear equation of state, and 
discussed energy and thermodynamic consistency issues. The above studies, however, did 
not clearly address the nature of the conversions between mechanical energy and internal 
energy, which is then addressed here in details. Section 2 provides the general theory. 
Section 3 applies the result to elucidating the thermodynamics of a Boussinesq fluid 
with a linear equation of state that has been widely used recently in numerical study of 
turbulent mixing, as well as in the context of horizontal convection. Section 4 discusses 
the case of the Boussinesq primitive equations currently used in numerical ocean models. 
Section 5 summarises and discusses some implications of the results. 
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2. Thermodynamically consistent and inconsistent 
Boussinesq/anelastic models 

2.1. Some specific examples of inconsistent Boussinesq models 

To help set up the context and the motivation for the present work, it is useful to provide 
some specific examples of energetically and thermodynamically inconsistent Boussinesq 
approximations, which have played and often continue to play a key role in the theoretical 
and numerical study of many fluid flows of interest in traditional and geophysical fluid 
dynamics. 

2.1.1. Boussinesq fluid with a linear equation of state 

The first Boussinesq model of interest is one that is meant to approximate fluids with 
an equation of state close to linear, retaining all physical processes apart from sound 
waves from the largest scales down to the molecular diffusive and dissipative scales. Its 
governing equations are: 

El + l VS p = -ao(p- P o) k + ^ {21) 

Dt po po 

V-v = 0, (2.2) 
DT 

— = «V 2 T, (2.3) 

p = Po [l- a(T - T )} , (2.4) 

where v = (u,v,w) is the three-dimensional velocity field, 5P = P — Pq is the pressure 
anomaly defined relative to the reference Boussinesq pressure Pq — —pQg Q z, p is the den- 
sity, T is the temperature, po and To are reference constant density and temperature, k 
is the unit normal vector pointing upwards in the direction opposite to gravity, and go is 
a nominal value of the acceleration of gravity. Although such a model is neither energet- 
ically nor thermodynamically consistent, as made clear in this paper, it has nevertheless 
been extensi vely used in recent theoreti c al discussion of the ener getic s of horizontal con- 
vectio n, e.g. JPaparella fc Young ( 2002 ). Wang fc Huang ( 20051) and Winters fc Young! 



as wel l as in discussing the energetics of turbulent mixing in stratified fluids by 



Winters fc al (1995). Moreover, such a model also forms the basis for numerous direct 



numerical simulations of stratified turbulence, in the sense that in such studies, both 
k and v are usually interpreted as representing the molecular values of diffusivity and 
viscosity respectively. From a thermodynamic viewpoint, the physical meaning of the 
temperature T is ambiguous, since it is conserved for adiabatic motions as if it were 
potential temperature, while also being homogenised by molecular diffusion as if it were 
in-situ temperature. Regarding pressure, it can be written as the sum of a purely hydro- 
static component plus a perturbation P' that can be regarded as the Lagrange multiplier 
associated with the incompressibility condition V • v = 0, and hence differs from the 
thermodynamic pressure. 

2.1.2. Coarse-grained Boussinesq model with parameterised turbulent fluxes 

A second important class of Boussinesq model is that associated with the kind of 
hydrostatic primitive equations model that has formed the basis for most numerical ocean 
models curren tly in used for climate change studies. An early formulation, following that 
introduced bv iBrvan" ( 1969h . is the following: 



^ + /k x u + iv^P = A ff V 2 ft u + A^, (2.5) 
Dt p oz z 
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V h • u + — = 0, (2.7) 

oz 

D6 , d 2 6 . . 

— = K H V 2 h e + K v -^, (2.8) 

-=A^5 + ^ — , (2.9) 

p = p(9,S,P), (2.10) 

where u = (it, u) is the horizontal velocity field, w is the vertical velocity, V/,. is the 
horizontal nabla operator, Ah and Ay represents horizontal and vertical turbulent eddy 
viscosities, while Kh and Ky represent horizontal and vertical turbulent eddy diffu- 
sivities. Here, the full nonlinear equation of state for seawater is used, and is usually 
formulated in terms of salinity S, potential temperature 9 and pressure P. Regarding the 
latter, it is worth pointing out that there is currently much debate about whether the 
Boussinesq pressure Pn = —pnqnz or the full hydros t atic p ressure P should be used in the 
equation of state, see Shchepetkin fc McW illiams | ( 201ll ) for a discussion of this point. 



Note that turbulent closure schemes in modern versions of primitive equa tions models are 



in general significantly more sophisticated than in the above model, see iGriffies I (|2004f> 
for a detailed discussion of current numerical ocean model formulations. Until very re- 
cently, it was generally thought impossible for Boussinesq primitive equations models 
such as the one above to admit a closed energy budget when using a nonlinear equation 
of state; as a result, numerical and empirical considerations must have been key in con- 
trolling the energetics of current numerical implementations of numerical ocean models 
and hence their numerical stability. 

2.2. Dynamical/ thermodynamic coupling in Boussinesq/ 'anelastic models 

The first step toward constructing an energetically and thermodynamically consistent 
Boussinesq/anelastic approximation for an arbitrary nonlinear equation of state is to 
clarify the nature of the coupling between the dynamics and thermodynamics in a Boussi- 
nesq/anelastic fluid. In a real fluid, such a coupling is achieved via the work of expan- 
sion/contraction which occurs through the term PDv/Dt, where P is the total pressure, 
and v the specific volume. To fix ideas, we examine the issue in the c ontext of the fol- 
lowi ng set of equati ons, which is based on that previously considered bv llngersoll (2005) 
and lPauluisI 1)20081 ) : 



^ + V.p)=-(^W + iv.S, (2.11) 
Dt V Po J \ Po J Po 



Drj 
151 





V • 


(pov) = 0, 


(2.12) 


= V = 


q _ 
T 


— V • {poF n ) + fjirr, 

Po 


(2.13) 


DS 

151 


= s-- 


= -— V-(poFs), 
Po 


(2.14) 




p = 


p(S,r},P ). 


(2.15) 



The above system of equations is sufficiently general that it also includes the usual 
Boussinesq approximation po = constant as a particular case, as is easily verified. For 
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that reason, Eqs. (|2. 11112.15]) will be hereafter referred to as the Boussinesq-Anelastic 
system, or BA system for short, where rj is the specific entropy, S is salinity, go is a 
constant acceleration of gravity, Z is the geopotential height, so that goZ = $ is the 
geopotential, and S is the stress tensor. The terms 77 = q/T and S are used as short-hand 
to denote the diabatic effects affecting 77 and S. These are further specified in terms 
of an entropy flux P„ and irreversible entropy production term rj irr > in Eq. (I2.13[) , 
which is the expression of the second law of thermodynamics, as well as in terms of a salt 
flux in eq. (|2.14[) . so that salt is assumed to be a conservative quantity. It is important 
to remark at this stage that the definition of buoyancy b — —go(p — po)/po is the one 
that is the most commonly encountere d in the l i teratu re. T his diffe r s, ho wever, from 
the form b = go(v — Vo)/vo assumed by Pauluis I ( 2008h a nd Young I 1 2010). As sho wn 
in this paper, it turns out that the form used bv IPauluis" ( 20081 ) and Young I ( 2010h is 
the one that should be used in practice, and the one that is the most energetically and 

thermodynamically consistent. 

Before proceeding, it seems important to point out that although IPauluis ] (l2008l) 's 
derivation of the above BA system initially assumes the reference pressure Pq to satisfy 
the classical hydrostatic balance p ~ 1 dPo/dz = —godZ/dz, this is no longer the case in 
the final energy conserving form of the equations. Indeed, in absence of fluid motion, 
setting P and p to zero in Eq. (|2.11l) shows that Pq and po must actually be solution of: 



d_ fPo 

dz \po 



■H7 — = -90 



8Z_ 

dz ' 



(2.16) 



which can be immediately integrated as: Po = P a — pogZ, where P a is the assumed 
constant and spatially uniform atmospheric pressure at Z = 0. This form is identical 
to the Boussinesq reference pressure, except for a non-constant po- This in turn implies 
that: 



D Po DPo 



Dt po 



Dt 



P0P0 



Dvo 
Dt 



v • (p v) = -pogoj^j- 



DZ DP, 



Oh 



Dt 



(2.17) 



where Poh is the reference pressure in hydrostatic equilibrium with po, which shows that 
the relative error in Po as compared to Poh, i-e., SPo/Pq = O(dpo/po), and is therefore 
small by assumption. Therefore, even though the anelastic approxim ation attempts to 
relax the incompressibility constraint of the Boussinesq approximation. IPauluis ( 20081) 's 
energy conserving form does so at the price of distorting the hydrostatic modes of mo- 
tions, so that it would be of interest to study the dynamical consequ ences of the BA 
system in more details, for instance by following iDavies &: et al. (2003)'s appro ach. The 
issue o f whether the distortion of the hydrostatic modes, which is not discussed in Pauluis I 
( 20081 ). could invalidate the approach is beyond the scope of this paper, which is primarily 
concerned with thermodynamic and energetic consistency issues. 

In order to discuss the energetics of the BA system, our starting point is the evolution 
equation for the kinetic energy, obtained in the usual way by multiplying Eq. (I2.11[) by 
/3qv, which after some manipulation can be written as follows: 



Po 



D_ v 2 
DtY 



V • (<JPv - poFfe e ) = p b- 



DZ 
Dt 



Po£a", 



(2.18) 



where the work against the stress tensor has been written as v • VS = V • (v • S) — po^K — 
V • {poFke) — Po£k as the difference between the divergence of a flux term minus the 
viscous dissipation term. For instance, in the special case where Pq 1 ^ ■ S = ^V 2 v, then 
F fce = ;;Vv 2 /2 and Sk = ^(I I Vm!! 2 + || V£p + || Vw|| 2 ) . Over the past decad e , a number 
of studies bv llngersoll I (|2005l) . IVallis I (|2006l ). IPauluis! (|2008l ). lYoung I (|2010f ). iNvcander " 
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(|2010h . and lMcIntvrel f|2010h all followed a similar approach in seeking to establish the 



energetic consistency of the above BA system, with the studies differing only in minor 
details and in the particular set of equations considered. Specifically, all these studies 
approached the problem by introducing the following function: 

,z 

h t {S,r l ,Z) = hl(S,r))- b(S,n,Z')dZ' (2.19) 
Jz 

obtained by vertically integrating the buoyancy b regarded as a function of Z and of the 
adiabatically conserved variables, so that: 

^ = - h ^ + C sS + Cr,V, (2-20) 

where 

dhl [ z dh x ,„, „ dhl ( z dh x 



Cs = i-/ 4 ^ c ' = w-/ a w dZ ' (2 - 21) 

Then, summing Eqs. ([2"TTgf and (pT2"0]) yields 

Po d1 (t + h *} + v ' (<5Pv ~ PoFke) = Po ( Cs ^ + C ^) - p° £k ■ ( 2 - 22 ) 



■R 



Because the right-hand side TZ of Eq. (I2.22j) obviously vanishes for adiabatic and inviscid 
motions, it may be concluded that the quantity po(v 2 /2 + h}) should be regarded as the 
relevant energy quantity being conserved in absence of diabatic and viscous effects. Note 
that the energy quantity thus constructed is no n-unique, since i t involves the arbitrary 



constant of integration hl(S,n) as discussed by iPauluis I (|2008j). In the context of the 



Boussinesq approximation, both I Young I (l2010h and lNvcander l (|2010l) chose Z = hi = 



and referred to h$ as the "dynamic enthalpy" and "effective potential energy" respec- 
tively. 

The above proof is incomplete, however, because establishing the energy consistency of 
the BA system actually requires: (1) Identifying the quantity playing the role of the total 
energy; (2) Demonstrating that the total energy thus identified is a conservative quantity 
for both diabatic and adiabatic motions. By contrast, all what the above arguments 
manage to establish is that the "adiabatic" energy v + h* is a conservative quantity 
for adiabatic motions only. When diabatic effects are retained, this adiabatic quantity is 
obviously created or destroyed owing to the r.h.s. of Eq. (|2.22l) becoming non-zero, and 
hence no longer conservative. In that case, the principle of energy conservation calls for 
the existence of an additional energy quantity h* obeying an equation of the type: 

pa^f + V -(poFh*) = -K, (2.23) 

involving some energy flux poF^ to be determined, as upon summation with Eq. (|2.22j) . 
the following conservation equation is obtained 

Po §t ( y + ht + h *) + v ' ((5Pv ~ poFke + po¥ * h) = °' (2 - 24) 

which now states that the energy quantity v 2 + h* + h* is conservative for both diabatic 
and adiabatic motions, and hence the natural candidate to play the role of the total 
energy in the BA system. 

Physically, the above approach casts the discussion of energetics in terms of the in- 
teractions between the adiabatic energy v 2 /2 + h^ and the diabatic energy h*. While 
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legitimate and s ometimes useful (as attested by the papers by iNvcander I (|201flh and 
iMcIntvrel ( 2010h for instance) , this approach cannot be regarded as entirely satisfactory, 
since: 1) the adiabatic and diabatic energies are fundamentally ill defined quantities, 
given that equally valid alternatives can be obtained by adding and subtracting any ar- 
bitrary function of S and rj from the former and to the later respectively; 2) the approach 
does not appear to be naturally linked to the classical description in terms of interactions 
between kinetic energy, gravitational potential energy and internal energy For instance, 
can we regard the sum h* + h* as equal to the sum of gravitational potential energy and 
internal energy, as for actual fully compressible fluids? How do energy conversions in the 
BA system compare with that of a fully compressible fluid? 

To address the above questions, a fundamentally different approach is clearly needed. 
If we are to understand the coupling between mechanical energy and internal energy in 
the BA system, it seems necessary to first identify what is the relevant expression for the 
mechanical energy, and then look for its evolution equation in order to determine how 
it couples with the internal energy. The most natural definition for the later appears to 
be: poE m — p [v 2 /2 — b(Z — Z )] = p w 2 jl + (p — p a )g {Z — Z ) (for some reference 
gcopotential height Z ), for which an evolution equation can be obtained from the kinetic 
energy equation (I2.18[) as follows: 



D 



— -b{Z- Z ) 



V • (5Pv - p Q F ke ) = -p (Z 



~Dt ~ p 



or equivalently, in conservative form: 



dt 



Pa 



5P 
Po 



PoFke 



-p (Z - Z ) 



Db 
Dt 



Pq£k- 



(2.25) 



(2.26) 



For comparison, the corresponding evolution equation for the mechanical energy in a 
fully compressible fluid takes the form: 

,2 



D 



V 

T 



g (Z - Z ) 



V ■ (Pv - pF ke ) = pP 



.Dv 
Dt 



pe K - 



(2.27) 



In a compressible fluid, mechanical energy and internal energy are coupled in mainly two 
ways: (i) via a reversible conversion between internal energy and kinetic energy accom- 
plished by the work of expansion/contraction PDv / Dt; (ii) the irreversible dissipation of 
kinetic energy into heat by viscous processes ek, as seen in the r.h.s. of Eq. (|2.27j) . Since 
viscous dissipation is obviously present in both Eqs. (|2.25[) and (|2.27|) . the only question 
that needs addressing is whether the following equivalence can be established? 



s Db Dv 
p Q (Z-Z Q )— O pP- 



? (2.28) 
Dt r Dt y J 

To see that this is indeed the case, we expand the pressure and specific volume as follows: 

J_ _ P - Po 
Po 



P = P„ 



Po9oZ 



Pi 



(2.29) 



which uses the underlying assumption in the Boussinesq and anelastic approximations 
that (p — po)/po *C 1. This implies in turn that at leading order, the work of expan- 
sion/contraction can be approximated by: 



PP 



Dv 
~Dt 



Po{P a ~ Po9oZ) 



D 
Dt 



1 

Po 



Po 



Pi 



= -Po(Z - Zo)^ + Po(Pa - pagoZ)^-. 



(2.30) 
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Eq. (|2.30p shows that the equivalence Eq. (|2.28[) is exact in the Boussinesq case po = 
constant, but only approximate in the anelastic case because of the term proportional 
to Dvo/Dt. This suggests, therefore, that the term SWba = {Z — Z$)db is the natural 
counterpart of the compressible work SW = —Pdv for a compressible fluid. The possi- 
bility to identify SWba with SW is a key result of this paper, for it points to a natural 
way to construct the whole range of known thermodynamic potentials for a BA fluid, 
as shown in the next subsection 12.31 To be fair, the identificatio n of SWba with SW i s 



not entirely new, since it wa s pointed out in earlier studies such as IWinters fc al I (|1995f ) 



Wang fe Huang 7 ! (|2005l) and lNvcander fc al. I (|2007l) for the kind of Boussinesq fluid with 



a linear equation of state discussed in Section [3] As far as we are aware, however, its 
consequences for constructing thermodynamically consistent Boussinesq and anelastic 
approximations do not appear to have been realized until now. Physically, Eq. (|2.28|) can 
be interpreted as expressing the fact that in a BA fluid, the compressible work manifests 
itself through the apparent changes in gravitational potential energy due to the apparent 
changes in the mass element dm = pdV , which is possible since a BA fluid conserves the 
reference mass element dmo = podV rather than dm. 

2.3. Thermodynamics of Boussinesq and Anelastic binary fluids 

Having established in the previous section that the energetics of both the Boussinesq and 
anelastic equations possess a term that can be identified as a conversion between internal 
energy and mechanical energy, it then becomes relatively straightforward to show that 
this can be used to construct energetically and thermodynamically consistent BA system 
of equations. To that end, let us first recall that a basic tenet of classical equilibrium 
thermodynamics is that the specific internal energy e = e(rj, S, v) of any fluid in local 
thermodynamic equilibrium (LTE) can be regarded as a function of state whose value is 
independent of the thermodynamic path followed, where rj is the specific entropy, S is 
salinity (the argument would als o work if S was the total mixing ratio, as for the case of 
moist air considered bv lPauluisI (120081) ). and v is the specific volume. As is well known, 



this implies that the reversible work transfer SW = —Pdv and generalised heat transfer 
5Q = Tdn + pdS entering the total differential of e, viz., 

de = 5Q + SW = Tdn + pdS - Pdv, (2.31) 

cannot be specified independently of each other. In classical equilibrium thermodynamics, 
this interdependence is imposed by the so-called Maxwell relationships, which simply 
express the result that the cross-derivatives with respect to two different variables must 
be equal for twice continously differentiable functions. In the present case, the assumption 
that e is a function of the thermodynamic state only implies that: 

de de de , 

~ 'drj' »=dS> P = -Tv> (2 ' 32) 
which in turn implies the following three Maxwell relationships: 

dT _ d 2 e dP dT _ d 2 e dp dp _ d 2 e dP 
dv dr]dv dn' dS drjdS drj' dv dSdv dS' 

The above remarks suggest that the simplest way to ensure that the BA system is en- 
ergetically and thermodynamically consistent is to ensure that the approximation to the 
generalised heat transfer, denoted here by SQba — Tb a di]b a + PbadSba, is similarly linked 
via relevant Maxwell relationships to the approximation to the work transfer identified 
previously, viz., 

SWba = {Z- Z )db, (2.34) 
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where the subscript ba was used to denote the approximation to the generalised 'heat 
variables and functions' T, 77, \i and S. If so, this would in turn allows one to regard the 
following expression: 

de ba = T ba dr] ba + /J, ba dS ba + (Z - Z )db, (2.35) 

as the natural counterpart of the fundamental relation of thermodynamics Eq. (I2.31[) . 
and hence as the total differential of the relevant approximation to the internal energy of 
the 'BA fluid', for which the natural variables are entropy 77, salinity S, and buoyancy b. 
In practice, however, it is often more convenient to work with pressure P as a dependent 
variable rather than specific volume, as well as with temperature T rather than entropy 
r), which motivates the introduction of additional thermodyn amic potentials constructed 
from Eq. (12.311) by means of the Legendre transform, e.g., lAlbertv I (|200lh . The most 



common thermodynamic potentials that are also the most relevant for the present work 
are the specific enthalpy h = e + pv and the Gibbs free energy g = e + pv — Tr], whose 
natural dependent variables are (r/,S,P) and (T,S,P) respectively. From Eq. (|2.35j) . it 
is easy to convince oneself that the corresponding approximations to h and g are given 
by: h ba = e ba - b(Z - Z ) and g ba = e ba - b(Z - Z ) - T ba T] ba , with the following total 
differentials: 

dh ba = d [e ba - b(Z - Z )] = T ba dr) ba + (J, ba dS ba - bdZ (2.36) 

dg ba = d [e ba - b(Z - Z ) - T ba i] ba ] = -rj ba dT ba + ^ ba dS ba - bdZ. (2.37) 

Until now, the above considerations have remained rather formal, and while they indi- 
cate that it is in principle possible to construct the full range of known thermodynamic 
potentials for a BA fluid, they have not addressed the issue of how such potentials might 
be constructed in practice. To simplify notations, we will usually refrain from using the 
suffix ba for the dependent variables in the following. We should keep in mind, however, 
that all thermodynamic variables used in the context of the Boussinesq/anelastic ap- 
proximations are fundamentally different, even if possibly only very slightly in numerical 
values, from their actual counterparts. 

In the p r evious approaches by llngersoll I ( 2005 ). Pauluis I ( 20081 ). lYoung I f|2010h and 



Nvcander 

I (l2010h . authors have generally assumed the buoyancy b = b(rj, S, Z) to be 



given as a function of the dependent variables considered, typically entropy, salinity and 
geopotential height Z. From Eq. (I2.36[) , this can be integrated with respect to Z, viz., 

h ba (r),S,Z) = h ba (r),S,Z )- [ b(r),S,Z')dZ', (2.38) 

for some reference geopotential height Zq, traditionally taken at the ocean surface in 
the oceanic case, but this only provides an expression for the specific enthalpy up to 
an indeterminate function of entropy and salinity, which cannot be specified without 
additional thermodynamic information about the fluid. 

The above problem arises because the knowledge of density alone is insufficient to 
specify all possible thermodynamic quantities; the heat capacity, for instance, cannot 
be inferred from density. On the other hand, many thermodynamic potentials have the 
property of e ncapsul a ting a ll known thermodynamic knowledge about a given fluid, as 
discussed by ICallen~l (|l985l ). In the following, we show how to obtain the relevant ap- 
proximations to the Boussinesq/anelastic system directly from two such thermodynamic 
potentials, the Gibbs function and the enthalpy. 
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2.3.1. Deriving BA thermodynamic potentials from exact thermodynamic potentials 

Given the similarity of the approximate and exact differentials for the internal energy, 
enthalpy, and Gibbs free energy, it is natural to ask whether it is possible to compute 
the approximate thermodynamic potentials from their exact counterparts, rather than by 
integrating the total differential of the approximate thermodynamic potentials. That this 
is indeed the case is shown first for the Gibbs free energy, which is the thermodynamic 
potential currently chosen to express all the thermod ynamic prop erties for seawater in 
th e latest standa rd recently adopted by UNESCO, see IOC (2010), as advocated earlier 



bv lFeistell (j2003h for instance. To that end, let us introduce the following function: 



9ba = 9ba{T, S, Z) = g(T, S, P oh {Z)) + g (Z - Z ) =g + g (Z - Z ), (2.39) 

and verify that its total differential agrees with the above differential relationship for 
the approximate Gibbs function, where g = g(T, S, P) is the exact Gibbs free energy 
expressed in terms of its natural variables T, S and P. Note that the approximate pressure 
used in Eq. (|2.39|) is the pressure Poh in hydrostatic balance with p , not the anelastic 
reference pressure P = P a — PogoZ. From the result that dg — —ndT + pdS + vdP, it 
follows that the total differential of gba is given by: 

dg ba = -fjdT + fidS + g (l - Po v)dZ, (2.40) 

where we used dPoh = —poQodZ, and where the tilde quantities are estimated for the 
reference pressure P h(Z) instead of the actual pressure, so that fj = r](T, S, P h(Z)), 
p = (j,(T, S, Poh(Z)) and v = v(T, S, P h(Z)). Eq. ([2~4TJ1) is obviously consistent with Eq. 
(|2.37[) provided that r/ ba , p ba and b be defined by: 

Vba = f-ba = P, b= -,9o(l - POV) = ~go ( ^ f° ) ■ (2-41) 



p 

The key result here is that while r]b a and fiba are found to be identical to their natural 
tilded expressions, the density p initially used to defined the buoyancy b = — go(p~ Po)/ Po 
is on the other hand found to differ from its natural tilded expression p, as Eq. (|2.41|) 
shows that: 

b=- 9 o(^A)=-9o( P -^), (2.42) 



which implies that: 



Po 

P = Po(l + ^y 2 ^) = Po (2 - Pov) . (2.43) 
Interestingly, the definition of the buoyan c y in E q. (12.41)) is identical to the one recently 



proposed by IPauluis I (|2008l ) and lYoung I (|2010l) . but differs from the definition usually 



adopted in most numerical ocean general circulation models as far as we know. 

As seen above, the knowledge of gb a is sufficient by itself to determine all possible 
thermodynamic potentials. Thus, the expressions for the specific internal energy and 
enthalpy are given by: 

e ha (S, T, Z) = g ba + b(Z - Z ) + T baVba =g+ {g + b)(Z - Z ) + Tfj, (2.44) 

h ba (S, T, Z) = e ba - b(Z - Z ) =g + g (Z - Z ) + Try, (2.45) 

which provide expressions for the specific internal energy and enthalpy in terms of the 
Gibbs function natural variables T, S, and Z, rather than in terms of e and /i's natural 
variables. It is easily verified that the total differential for e ba and h ba are given by: 

de ba = Tdf) + /idS + (Z- Z )db, (2.46) 
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dh ba = Tdfj + /idS* - bdZ, (2.47) 

which can be checked to be in agreement with Eqs. (I2.35P and (|2.36[) . thus completing 
the proof. 

While the use of the Gibbs function is the most natural approach when the thermo- 
dynamic properti es are tabula ted as functions of T, S, and P, as is currently the case 
for seawater, e.g., IOC ( 201dl ). it is sometimes more convenient, depending on the par- 



ticular situation considered, to work with different dependent variables and hence with 
a different "master" thermodynamic potential. For a general discussion of the d i fferen t 
ways to compile thermodynamic data, the interested reader is referred to C alien I ( 19851 ). 



The specific enthalpy h — h(rj, S, P), owing to the particularly important role it plays 
for describing and understanding the energetics of stratified fluids, warrants a separate 
discussion and is therefore discussed next. To that end, let us introduce the function h ba : 



T-ba 



h{t], S, P oh (Z)) + g a (Z -Z )=h + g (Z-Z ), 



and verify that it is the relevant Boussinesq/ anelastic approximation to the specific en- 
thalpy. This is done by taking the total differential of h ba , viz., 

dh ba = fdr) + jldS + g (l - p v)dZ = fdrj + fidS - bdZ, (2.49) 

and checking that it agrees with Eq. (|2.36l) . QED. The corresponding expression for the 
internal energy e ba = h ba + b{Z — Zq) can be written as: 

e ba = h+(g + b)(Z-Z ). (2.50) 

Its total differential is: 

de ba = fdT)+p,dS-p g vdZ+(b+g )dZ+(Z-Zo)db = fdr)+p,dS+(Z-Z )db, (2.51) 

which is again consistent with Eq. (|2.35j) . 

As a final remark, let us mention that while the above derivations demonstrate that 
the BA thermodynamic potentials can be constructed from the knowledge of the exact 
Gibbs function or enthalpy, integration of the Maxwell relationships will in general be 
needed to construct the BA thermodynamic potentials when an idealised equation of 
state is assumed. Sections [3] and |4] illustrate in details the two different possible cases. 

2.4. Consequence for energy conservation and "heat" -related quantities 

The above results make it possible to write down the evolution equations for the specific 
internal energy and enthalpy by combining Eqs. (|2.46[) and (|2.47[) with the evolution 
equations Eqs. (|2.13[) and (|2.14j) for 77 and S respectively, leading to: 

Po^ = p [q + fiS]+Po(Z-Zo)^ t , (2.52) 



Dh ba 



fiS 



DZ 

p b—. (2.53) 



Interestingly, taking the difference between these two equations yields the evolution equa- 
tion for the gravitational potential energy b(Z — Zq): 

D fur 7 w D(h ba - e ba ) ,„,.> 

Po Dt ^ ~ ^ = P ° Dt ' ^ ' 

which shows that for a BA fluid, the gravitational potential energy can be regarded as a 
thermodynamic property of the fluid, since it is entirely determined from the knowledge 
of S, T, and Z. There are therefore two main ways to describe the coupling between 
dynamics and thermodynamics in a BA fluid, the first one focusing on the coupling 
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between internal energy and the total mechanical energy (the sum of kinetic energy and 
gravitational potential energy), and the other between the kinetic energy and the total 
potential energy (the sum of internal energy and gravitational potential energy, which is 
equal to the enthalpy here). By combining either one of Eqs. (|2.52p or (12.53)) with the 
mechanical energy equation Eq. (|2.25|) , the two following equivalent evolution equations 
for the total energy E t = v 2 /2 + e ba — b(Z — Z ) = v 2 + hb a are obtained: 



Po 



D 

Hi 



b(Z - Z Q ) + e ba 



V • [&Pv - poFke 



po 



q + flS - e K 



(2.55) 



D 

P0 Wt 



V 7 

-jj" + h ba 



V • [SPv - poFfee] = Po q + p.S - e K 



(2.56) 



Now, in order for the BA system to be energetically consistent, it remains to verify that 
Eqs. (|2.55p and (I2.56|) are consistent with the principle of energy conservation. This is 
easily shown to be the case only if the right-hand side of Eq. (|2.55l) is the divergence of 
some energy flux F q , i.e., 



p (q + fiS- e K ) = -V • (p F q ). 



(2.57) 



As we show now, this imposes a constraint on the form of the entropy flux F,, and irre- 
versible entropy production r\ irr in Eq. (|2.13[) . Assuming salt to be conservative quantity, 
and hence such that poS = — V ■ (poFs) for some salt flux Fg, implies for the evolution 
equation of specific entropy: 



Po- 



Di%a _ Poq_ _ PqSk + pV ' (PoFg) - V ■ {ppF q ) 
Dt ~ T T 



Pq£k 
T 



p Q (F q ~ pFj 
T 



Po 



■v,i -,..v a 



(2.58) 



The latter expression establishes that the entropy flux F, ( and irreversible entropy pro- 
duction r]i rr in Eq. (|2.13p must be related to the salt flux Fg and enthalpy /internal 
energy flux F g by: 



F — 



F q - pFt 
T 



(2.59) 



f, irr = ^ + F q -v(±) -Pg-vf^ 



(2.60) 



In non-equilibrium thermodynamics, this is usually the point at which the second law of 
thermodynamics is th en invoked to furthe r const rain F q and Fg, in order to guarantee 
that fjirr ^ 0, e.g., see lde Groot &: Mazur ] (Il962h . Going back to the evolution equation 
for internal energy and enthalpy, note that Eq. (|2.57|) implies: 



Po-Jjjr = ~ V ' (Po F q) + Po^k + p (Z - Z )^, 



(2.61) 



Po- 



Dh ba 



-V • (p F q ) + pqSk - Pob 



DZ 



(2.62) 



Dt ' vru H ' ' r "" Dt 

so that F g appears as the diffusive flux of internal energy or enthalpy. Both equations 
can be regarded as equivalent forms of the first law of thermodynamics. 
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2.5. Alternative forms of the first law of thermodynamics 

We now investigate the consequences of the above results for the determination of various 
temperature variables. Let us first note that from the differential of the Gibbs function: 

d<?ba = — fj&T + fidS — bdZ, 

the Maxwell relationships provide the following partial derivatives for the entropy: 

dfj dfi dfj db dv 

^s^'dr' ~dz = df = Po9o df = 

by defining a = podv/dT as the relevant definition of the thermal expansion coefficient. 
Using the well known result that drj/dT = c p dT/T, this implies that we can write: 

dr? = C ^dT - ^dS + 90 &dZ (2.63) 

As a result, it follows that the temperature equation can be written as: 
DT T Drj T dp, DS ag G T DZ 



Dt c p Dt c p dT Dt c p Dt 

In absence of salinity, this equation takes the simpler form: 

DT _ V • (KCpVT) e K ag T DZ 
Dt Dt ■ 



(2.64) 



(2.65) 



by assuming the diffusive heating to be given by the classical Fourier law V • (poF q ) — 
— V-(KpoCpVT). This shows that the evolution equation for in-situ temperature in general 
possesses: a) a term related to molecular diffusion, that is general not-conservative, i.e., 
it cannot be expressed as the divergence of a flux because c p is not constant; b) it usually 
incorporate the Joule heating due to viscous d issipation; c) it possesses a term related 
to change in pressure, which some authors, e.g. lPons fc Le Querel (|2005l . 120071 ). call the 
"piston effect" , and the resultant Boussinesq equations, the thermodynamic Boussinesq 
equations. 

In practice, the pressure effect can be accounted for by constructing evolution equa- 
tions for the potential temperature or conservative temperature respectively. In order to 
show how an equation for can be constructed, it is useful to construct the potential 
temperature 9 first. By definition, 9 is the temperature that a parcel with temperature 
T would have if lifted adiabatically to a reference level Z = 0. Potential temperature is 
thus defined as the implicit solution of the following equation: 

Vba(T, S, Z) = r)ba(9, S, 0). (2.66) 
Differentiating this expression yields: 

d^ a = |dT - ||dS + a go dZ = °jd9 - ||dS, (2.67) 

where we defined c r p = c p (9,S,0) and fi r = fl(9,S,0). In the following, we simplify 
notations by defining c r p = c p (9,0). This expression shows that it is possible to rewrite 
the evolution equation for entropy as follows: 



c^D9__dpm Dn ba 
~9~Dt ~ dS~Dt + Dt 
In absence of salinity, we can write: 

c r p D9 _ V • (ncpVT) | e K 



(2.68) 



9 Dt T T 



+ -^- (2-69) 
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This can be transformed in the following equation for 9, 



D9 



= V- 



Kc p evT\ 



Dt V c r p T J 

where the nonconservative production/destruction of 9 is given by: 



%rr — KCpS/T • V . 



8e K 



(2.70) 



(2.71) 



In the literature, 9 is often treated as a conservative variable, which consists in neglecting 
the nonconservative term # irr . Alternatively, one may remark that the equation for 9 can 
also be written as: 

~ r D9 Dh e , 9 79 s 

where hg{9) = hb{r], 0) is the enthalpy a parcel would have if move d adiabatical l y from 
Z to Z = 0. For this reason, hg was called "potential enthalpy" by iMcDougall I ()2003h . 
This allows one to defined a new temperature variable 0, also conserved for adiabatic 
motions, such that: dhg = c° d6, where c° is an arbitrarily defined specific heat capacity 
representative of c p at Z = 0. IMcDougall I (|2003l ) discusses a possible choice for c° in the 
oceanic context. As a result, it is possible to write the above equation as: 



DO 
Dt 



1 Dh 6 
$~Dt 



= V 



c°T 



where this time, the nonconservative production/destruction of & is given by: 

9s K 
T 



KC 

'~c\ 



(2.73) 



(2.74) 



McDougall I (|20Q3I ) provides convincing evidence that the nonconservative term <9j rr is 
significantly smaller than 9i rr for the present-day oceans. The relative smallness of Qi rr 
over 9i„ does not appear to be entirely generic though, and should therefore be checked 
on a case by case basis, as it may depend on the particular fluid considered and on the 
degree of turbulence present. It is beyond the scope to simplify the equations further, 
and to discuss when the evolution equation for O can be approximated by the above 
simple diffusive model using kV 2 T. 



3. Thermodynamically and energetically consistent model for a 
Boussinesq fluid with a linear equation of state 

In this section, we return to the case of a Boussinesq fluid with a linear equation of 
state, originally assumed to be governed by Eqs. (l2.HI2.4j) . with the aim of clarifying its 
thermodynamics, as well as to improve its formulation to make it thermodynamically 
and energetically consistent. 

3.1. Thermodynamic properties of a Boussinesq fluid with a linear equation of state 

We first seek to construct the relevant forms of the specific Gibbs function, enthalpy and 
internal energy, assuming a linear equation of state p — p^[l — a(T — To)] and constant 
specific heat capacity c p o, so that the expression for the buoyancy is: 



b = ago(T - 7n). 



(3.1) 
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We need to decide on whether to interpret T as the potential or in-situ temperature in the 
equation of state, as we saw above that any thermodynamically consistent formulation 
must distinguish between the two kinds of temperature. For completeness, the two cases 
are considered. We first discuss the case where T is the in-situ temperature, deferring the 
discussion of an equation of state linear in potential temperature to subsection 13.31 This 
choice amounts to assume that the isothermal compressibility, rather than the adiabatic 
compressibility, vanishes. From a thermodynamic viewpoint, this is arguably an awkward 
assumption, because Eq. ([l.ip suggest that this endows the fluid with a negative speed 
of sound, but since sound waves are filtered out, it is unclear whether this is a serious 
impediment. 

The interpretation of T as the in-situ temperature motivates looking for a description 
of the thermodynamic properties in terms of the specific Gibbs function g ba , since the 
natural variables of the latter are T and Z. As seen previously, the total differential of 
g ba is dg ba = -r) ba dT - bdZ, which implies: 

dgba dg ba 

~~dT~ = - Vba ' ~dZ = ^ (3 ' 2) 
The assumption of constant heat capacity provides the additional differential relation: 

rr d 2 gba ,„ „n 

The system of partial differential equations Eqs. (|3.2j) and (|3.3[) is sufficient to completely 
specify g ba , whose integration yields 

9ba = -c p0 [Tin (T/T ) - T] — ag (Z - Z )(T - T ), (3.4) 

up to some arbitrary constant of integration. From Eq. (|3.2p . the following expression for 
the specific entropy is obtained: 

Vba = = Cpo ln (Jfr) +&9o{Z - Z ). (3.5) 

This makes it possible to derive an exact expression for the potential temperature 9, 
which is the implicit solution of the equation r] ba (T,Z) = rj ba (9,0), i.e., c p oln(T/T0) + 
ag (Z - Zq) = Cp ln {9/T ) - ag Z a , yielding: 

tf = T exp(^l (3.6) 



CpO 

This makes it clear that in order for a model to be thermodynamically consistent, po- 
tential temperature is always different from the in-situ temperature. Now, we previ- 
ously established that the expressions for internal energy and enthalpy were given by: 
f ba = 9ba + b(Z - Zq) + Ti] ba and h ba = e ba - b(Z - Z ) = g ba + Tr/ ba , yielding therefore 
the following expressions: 

h ba — Cp a T + ag T (Z - Z ), (3.7) 
ago(Z - Z Q ) 



eba = c p o 



1 



T. (3.8) 



CpO 

These relations can be written in terms of other variables by using the expression for 9, as 
well as b. For instance, the enthalpy can be written in terms of the potential temperature 
and Z as follows: 

ha = c p0 9 exp ( - a9 ° Z | + ag T (Z - Z ). (3.9) 



CpO 
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^ _ agoZp ^ J ag Q 



c p0 [b + ag T 



T + —) (3.10) 



ago 



T = e^ P i-^\ (3.15) 



3.2. Improvement of the model energetic and thermodynamic consistency 

Having clarified the nature of the thermodynamics and thermodynamic potentials for a 
Boussinesq fluid with a linear equation, we now turn to the issue of writing down explic- 
itly how to transform the initial formulation into an energetically and thermodynamically 
consistent one. The main modification introduced is related to the heat equation. More- 
over, there is now a distinction between in-situ temperature and potential temperature. 
To summarise, the whole model is therefore given as follows: 

— + — V5P = bk + vV 2 v (3.11) 
Dt po 

b = ag (T-T ) (3.12) 

V-v = (3.13) 

} 

CpQ J 

Note that we can rewrite the temperature equation as: 

^ = V.(^VrW + ^ (3-16) 
Dt \T J Cp T v ' 

where 

6 irr = — tVT ■ V (1) = ^ = + ^0 (3.17) 

\T J c p o T oz c p0 dz c* 

Note that it is also possible to write the flux entirely in terms of 0: 

^VT = «V0-^k (3.18) 
T c pQ 

which was computed assuming Z(z) — z. Interestingly, these equations introduce the 
following length scale: L — c p o/(ago)- For typical values, c p o = 4.10 3 J.kg _1 .K~ 1 , a = 
10 _4 K~ 1 and go = 10m.s — 1 , which yields: L = 4.10 6 m. As L is huge compared to 
the typical length scales at which molecular diffusion is important, it follows that the 
nonconservative term is probably negligible. 

3.3. Equation of state linear in potential temperature 

This section revisits the above results by assuming that the equation of state is linear 
in rather than in T, and hence that the adiabatic compressibility vanishes, as is gen- 
erally implicitly assumed in traditional low Mach number asymptotics. In this case, the 
buoyancy becomes: 

b = a e go{0- o ), (3.19) 
with ag the isentropic thermal expansion, which Tailleux 1 ( 2010f) shows (his Eq. A. 11) is 



related to the classical thermal expansion by ag — aT/0 for a constant c p q. The use of 
as a dependent variable makes it natural to seek a description of the thermodynamic 
properties of the fluid from the specific enthalpy, since is closely related to the specific 
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entropy by d?7b a = c p od9/9, as seen previously. As a result, the total differential of specific 
enthalpy can be written as: 

dh ba = Tdr] ba - bdZ = ^^-d9 - bdZ, (3.20) 

9 

which yields the following system of two partial differential equations: 

dh ba c p0 T dh ba . . 

~qq = ~q ' -Qz- = -b=-a e g {8-9 () ). (3.21) 

The Maxwell relationship for such a system imposes the following integrability constraint: 

9 fc p0 T\ db 

dz{—) = -de = - ae90 > (3 - 22) 

which can be integrated to yield the following expression between T and 9: 

T =1 _ aegoZ^ ^ 
c p o 

using the fact that by definition, T — 9 at Z — 0. By using Eq. (|3.23p . it becomes 
straightforward to integrate the partial differential relations Eqs. p.21j) to eventually 
obtain 

a s g Z 

CpO 



' - {0 -0 )=[ ) c p0 T, (3.24) 

up to an arbitrary constant of integration. The expression for the internal energy is 
therefore given by: 

e&a = h ba + b(Z - Zq) = [c p0 - a e g Z }(9 - 9 ), (3.25) 

which implies that in the present model, 9 is in fact a proxy for the specific internal 
energy. It is useful to note here that because c p o is constant, the concepts of potential 
temperature 9 and conservative temperature O coincide. 

The evolution equation for temperature is unchanged, while the expression linking 
T and 9 is now given by Eq. Q3.23p . so that in summary, the thermodynamics is now 
described by the system: 

DO 9 9 9e K , 
— = -kV 2 T + — -, (3.26) 
Dt T c p qT 



T = 



(3.27) 



c P o 

As shown in this paper, the diffusive term can be written as the sum of the divergence of 
a diffusive flux, plus a nonconservative production/destruction term, which are explicitly 
given by: 

™ = v .(| VT ) +9 ,„ + H ( , 28) 

where 

\Tj [l-aegoZ/c p o} 2 dz c p0 \T J dz 

The nonconservative terms have been shown by iMcDougall ( 20031 ) to be very small 
compared to the diffusive term, so that can be neglected in practice. 
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4. Thermodynamically and energetically consistent Boussinesq 
primitive equation ocean models 

In this section, we return to the Boussinesq primitive equations Eqs. (j2.5H2.10j) that 
form the basis for most current numerical ocean general circulation models (save for the 
form of the turbulent parameterisations) , with the aim of showing how to modify them 
in order to make them energetically and thermodynamically consistent. 

4.1. Improved formulation 

It should be clear by now from the above results that two main modifications are needed to 
make the Boussinesq primitive equations (j2.5H2.10]) energetically and thermodynamically 
consistent, which are: 1) modification of the definition of buoyancy, 2) addition of the 
nonconservative terms in the thermodynamic equation Eq. (|2.8p . The resulting set of 
modified equations is as follows: 

Du 1 „ o d 2 u 

— + /kxu+— V h 5P = A H Vlu + Av — , (4.1) 

i«p JP^yz^ez 



po dz \ p J dz dz 

d't 

dz 



Dw 

V h -u+^=0, (4.3) 



DO d 2 B 

— = K H v 2 h e + K v — + e irr , (4.4) 

Vt oz z 

™=K B ViS + K 9 ™, (4.5) 

p = p(S,Q,P (Z)), (4.6) 
where 9 is now assumed to be the conservative temperature, on the basis o f lMcDougalll 



(|2003l )'s results suggesting that the nonconservative term Qi rr is significantly smaller for 
than for the potential temperature 9, resulting in a smaller error in the overall energy 
budget when this term is neglected. For completeness, we need to provide an explicit 
expression for Qi rr . To that end, it is first necessary to understand how to evaluate a 
number of thermodynamic properties for Boussinesq seawater, which is discussed next. 

4.2. Thermodynamics of Boussinesq seawater 

In contrast to the idealised Boussinesq model of the previous section, the aim of numerical 
ocean modelling is to achieve realistic simulations, and therefore to use an equation of 
state that is as realistic as possible. As mentioned earlier, the ther modynamic pro p erties 
of sea water are currently encapsulated in a Gibbs function, e.g., iFeistel I (|200ah :llOC' 



( 20101) . with natural variables in-situ temperature, salinity, and pressure. It is therefore 



natural to seek a determination of Boussinesq thermodynamics from the published Gibbs 
function. The Boussinesq Gibbs function was shown above to be given by: 

g ba (T, S, Z) = g(T, S, P (Z)) + g (Z - Z ) =g + g (Z - Z ). (4.7) 

As is well known, the specific entropy r]b a , the relative chemical potential /iba-, the buoy- 
ancy b (and specific volume), and the specific heat capacity c Pt b a can all be obtained from 
the first and second partial derivatives of g ba as follows: 

Vba (T,S,Z)=fj=-^;, (4.8) 



22 



R. Tailleux 

Vba{T,S,Z) = p,= ||, 



dgba 
dZ 



dl 
dZ 



9a = 9o(PoV 
_ T d 2 9 



(4.9) 
(4.10) 
(4.11) 



^pa — ^ - - Qj, 2 

where the formula for b was arrived at by using the result that dg/dP = v(S, T, P). As 
previously, tilded quantities refer to quantities estimated for the hydrostatic reference 
pressure Po(Z) rather than the full pressure, so that v = v(T, S, Pq(Z)) for instance. 

In order to arrive at an expression for the conservative temperature, we need the 
expressions for the specific enthalpy h ba and potential enthalpy he, which are respectively 
given by: 

h ba (T, S, Z) = g ba (T, S, Z)+Ti) = g + g (Z - Z ) + Trj, (4.12) 

h e (T, S, Z) = h ba (6, S, 0) = h r (9, 5), (4.13) 

where the potential temperature 9 is now the implicit solution olr]i ia {T 1 S, Z) = rj ba {9, S, 0). 
Using the above expressions for the specific entropy and heat capacity, as well as the 
Maxwell relationships for the Gibbs function, yields the following expression for the total 
differential of specific entropy: 



drjba = d?7 



CpdT 



dS 



db 



dZ, 



(4.14) 



T dT dT 

which in turn yields the following expression for the total differential of specific enthalpy: 



dhha = Tdf] + jldS — bdZ = T 



_ dT dfl 



dS- 



db_ 



dZ 



-fidS — bdZ 



di] 



CpdT - 



fL-T^)dS+(g &T-b)dZ, 



(4.15) 



where db/dT = p^g^dv/dT = go&, by defining the thermal expansion coefficient by 
a = podv/dT. Similarly, the haline contraction coefficient should be defined by /3 = 
—podv/dS. By evaluating the above expression at Z = 0, we obtain the following ex- 
pression for the total differential of potential enthalpy: 



dhe = cldO 



^ ) ds = c °de. 



(4.16) 



where c r p = c p {9, S 1 , 0) and fi r = fl(9, S, 0). As discussed by iMcDougall I (|2003l) . the con- 
servative temperature is defined such that d0 = dhg/cp, for some constant representative 
value of the specific heat capacity c°. The latter expression allows one to express d9 in 
terms of d6 as follows: 



c u 1 
d6 = JL d <3 - ± 



p, 



- 9 d £\ dS. 



(4.17) 



Now, by evaluating the differential of entropy Eq. (|4.14j) at Z = 0, as well as by using the 
implicit definition of potential temperature and Eq. (|4.17[) . it is possible to write down 
the differential of entropy equivalently in terms of 9 or O as follows: 



(4.18) 
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which in turn allows one to write the total differential of specific enthalpy in terms of 
the conservative temperature as follows: 



Tc° 

dhba = Tdf] + pdS - bdZ = -^-d@ 



Tfl r 



dS - bdZ. 



(4.19) 



This result is used next to provide an explicit expression for the nonconservative produc- 
tion of conservative temperature. 

4.3. Determination of the nonconservative term !rr 



As discussed in iTailleux I (|2010D . the determination of @i rr follows from the constraint 
of total energy conservation. In order to see this, let us first derive the evolution equa- 
tion for the kinetic energy by multiplying the momentum equations by v. After some 
manipulation, this can be put in the form: 

D u 2 DZ 

Pa ~blT + v ' (<5Pv ~ PoFke) = pah ~b~t ~ poSK ' (4 ' 20) 

where the diffusive flux of kinetic energy and viscous dissipation are given by: 

F fce =A H V h ^-+ A v ^-^r-k, (4.21) 
2 dz 2 



e K = A h \\\V h u\\ 



A, 



dz 



(4.22) 



Now, by inserting the evolution equations for O and S into Eq. (|4.19j) . it is possible to 
rewrite the evolution equation for the specific enthalpy as follows: 



po- 



Dh ba Tel DO 



Dt 



Tfj, r \ DS 



DZ 



Dt - Pob -Dt 



Tc° ■ DZ 

V • (poF/J - Poh irr + —^-poQirr ~ PV b ~f^ ') 



(4.23) 

where the diffusive flux F/, and the irreversible nonconservative term hi rr are given by: 



K H V h e + K v — k 
dz 



H 



//.,. 



Tc2 



Tp r 



K H V h S + K v —k 
oz 



(4.24) 



K v 



dS dp** dQ d ( Tc° p s 
dz dz dz dz \ 9 



(4-25) 

where p** = p — Tp r /9. Finally, an evolution equation for the total energy is obtained 
by Eqs. @jZDH and flgggp , 



Po 



D (y 
Dt 



( . Tc° . \ 

V ■ (SPV + p F h - p F ke ) = po h irr + —^-&irr ~ Po£ K , (4.26) 



which shows that the total energy u 2 /2 + h\> a is a conservative quantity only if the 
right-hand side of the later equation vanishes, which yields: 



@7>r — 



£k + hi 



(4.27) 



This completes the full specification of the modified Boussinesq primitive equations. 
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4.4. Errors in the energy budget of current OGCMs 

The above results imply that there are at least two main sources of error in the global 
energy budget of current OGCMs, resulting respectively from their use of the buoyancy 
bogcm = —go{p — Pq)/po instead of b = —go(p — Po)/p, as weu as from neglecting the 



irreversible production term &i rr in the temperature equation (see I Tailleux I (|2010f ) for 



relevant expressions when potential temperature is used). Thus, by replacing the buoy- 
ancy b by boccM in the kinetic energy equation, as well as by neglecting Oi rr in the 
enthalpy equation, implies for the evolution equation of total energy in current OGCMs 
to be of the form: 

D /v 2 \ 

Pa Dt [ Y + hba ) + v ' ^ Pv + PoF '' ~ PoFke ^ = poEirr > ( 4 - 28 ^ 

where the right-hand side p$E irr corresponds to the spurious production/destruction of 
total energy resulting from the inconsistent treatment of the buoyancy and temperature 
equation. From the above results, its explicit expression is given by: 

p Q E lrr = -p (hirr + £k) +Po(p ogcm -b)w « -p (ji irr + e K j - - P ° ^j p w, (4.29) 
since we have: 

b ogcm -b = -g ( tlfl) +90 ( Ll*) = go (2-I : -El) K} -( £Z^y . (4.30) 



Po J \ P J \ Po P / \ Po 

Note her e that the final form for the error term Ei rr differs significantly from that dis- 
cussed in lTailleuxl(|2010h . because in the latter study, Ei rr was estimated using the non- 



approximated thermodynamic potentials, resulting in additional spurious sinks/sources 
of energy. Interestingly, it is important to note that because the term (b ogcm — b) scales 
as the square of the small Boussinesq parameter (p — po)/po, it follows that the ener- 
getically and thermodynamically consistent Boussinesq formulation has the same order 
of accuracy as the non-consistent formulation. From that viewpoint, it can be said that 
of all the possible Boussinesq approximations that can be constructed with the same 
formal order of accuracy, one can be found that is energetically and thermodynamically 
consistent. 

As discussed in Tailleux ( 2010l ). errors in the global energy are likely to be of con- 



sequenc e only if they ar e associated with spurious forces in the momentum equations. 
Indeed, Tailleux ( 2010l ) estimated that the same error O(10~ 6 W.m -3 ) in the overall 



energy budget corresponds only to a spurious heat source/sink in the temperature equa- 
tion that is 0(7.5iT/(million years)), which is arguably utterly small, whereas it may 
correspond to spurious positive or negative acceleration in the momentum equations 
0(3m.s _1 /year) in the momentum equations, which is in contrast a significant number. 
On this basis, it appears important to correct the definition of the buoyancy in cur- 
rent OGCMs, but less vital to retain the nonconservative term Qi rr . From a practical 
computational viewpoint, modifying the definition of buoyancy in current OGCM im- 
plementations should be rather straightforward, as it amounts to change the integrand 
in the particular sub-routine estimating the hydrostatic pressure from the knowledge 
of density. By contrast, evaluating the nonconservative term Qi rr at each grid point at 
all time steps would be a major undertaking, since it requires the evaluation of many 
thermodynamic quantities such as p,, p r , c r p , 9 and T. This would be computationally 
prohibitive, given that the evaluation of density alone already amounts for a significant 
fraction of the total CPU time owing to its strongly nonlinear character. 
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5. Summary and discussion 

In this paper, we sought to clarify the nature of the conversions between mechanical 
energy and internal energy supported by the Boussinesq and anelastic approximations, 
in the general case of a binary fluid with an arbitrary nonlinear equation of state. A key 
result was to show that the energetics of such approximations possesses a term that can be 
identified as playing the role of the classical compressible work of expansion/contraction, 
which manifests itself as apparent changes in gravitational potential energy due changes 
in density (and hence of mass if the volume is exactly or approximately conserved, as is 
the case for such approximations). In contrast with a fully compressible fluid, however, 
the conversion between mechanical energy and internal energy is between GPE and IE 
rather than between KE and IE. By regarding this term as the Boussinesq/anelastic 
approximation to the compressible work 5W = —Pdv, it is possible to construct the 
relevant approximation to the "heat" SQ in a consistent way by ensuring satisfaction 
of Maxwell relationships, which upon integration eventually leads to the construction 
of consistent expressions for the specific internal energy, a s well as of the whole range 



of known thermodynamic potentials, illustrating iBannon I (|1996l )'s statement that the 
"thermodynamics is slaved to the dynamics" . 

The existence of well-defined thermodynamic potentials from which to derive physi- 
cally consistent expressions for the first law of thermodynamics appears to be sufficient 
to endow the Boussinesq and anelastic approximations considered in this paper with 
fully consistent energetics and thermodynamics, even when diabatic effects and a fully 
nonlinear equation of state for a binary fluid are retained. In the energetically consis- 
tent form of the Boussinesq and anelastic approximations, the sum of kinetic energy and 
enthalpy is a conservative quantity, and hence the natural total energy for the system. 
As a consequence, the gravitational potential energy can be regarded as the difference 
between enthalpy and internal energy, and hence as a purely thermodynamic property of 
the fluid. These ideas were illustrated by showing how two widely used but energetically 
inconsistent Boussinesq models could be modified to make them fully energetically and 
thermodynamically consistent. Interestingly, we find that the modifications required to 
ensure energetic and thermodynamic consistency do not alter the formal order of accu- 
racy of the approximations. In other words, in the space of all Boussinesq and anelastic 
approximations of a given order of accuracy, one exists that is fully consistent energet- 
ically and thermodynamically. It was also showed how to construct explicitly the full 
range of thermodynamic potentials for Boussinesq and anelastic models, either by inte- 
grating Maxwell relationships in the context of idealised models with idealised equations 
of state and heat capacities, or by approximating the exact thermodynamic potentials 
when those are known. A direct application of our results is to suggest that current nu- 
merical ocean circulation models possess a potentially significant source of error in their 
momentum equations owing to their use of an incorrect definition of buoyancy, which 
could in principle be simply corrected by using the correct definition. Interesti ngly, the 



impro ved d efinition of bu oyancy is one that was recently proposed earlier by iPauluis 



(|2f)f)8h and [Young] (|2010h . Dut it does not appear to have been realized until now that 



such a modification was needed to improve the energetic consistency of current OGCMs. 
An important implication of our results is to support the earlier suggestion bv lTailleux 



(2009) that the Boussinesq (and hence anelastic) approximations can support large con- 



versions between mechanical energy and internal energy, and therefore compres sibility ef- 
fects s ignificantly larger than previously assumed. Specifically, the point made in lTailleux 



(2009) is that in the context of turbulent stratified mixing, the apparent irreversible con- 



version of available GPE into background GPE should not be interpreted as a mechanical 
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to mechanical energy conversion, as proposed by IWinters fc al (1995), but as an irre- 



versib le conversion of AGPE into IE. Note, indeed, that on the one hand, IWinters fc al 
(1995) interpret the following energy conversion term: 



r.i'l III I IUI V 



g ZKV 2 pdV = ng [(p) bo tto 



(P) 



topi 



(5.1) 



as a (laminar) conversion of internal energy into background gravitational potential en- 
ergy, but on the other hand, interpret the following energy conversion term: 

gnWS/pr" 



W r . 



turbulent 



gz r nXl z pdV - W r , 



■ aminar 



dp r /dz 



-dV - W, 



r.la in i an r • 



(5.2) 



as the irreversible conversion of AGPE into background GPE r , where z r — z r (x, t) is the 
parcel's position in Lorenz's reference state, and hence as a mechanical to mechanical 
energy conversion, on the grounds that W ri turbulent shows up in the evolution equations 
for AGPE and GPE r with opposite signs. Yet, both W r> i am i nar and W r turbulent are seen 
to involve terms of the form gznV 2 p, which the present paper argues is the one playing 
the role of the compressible wo rk of expansion/ contraction in the Boussinesq/anelastic 
approximations. For this reason. iTailleux I (l2009h argue that W ritU rbuient actually refers to 
two different types of energy conversions in the AG PE and GPE r evolution equations, 

~ (l2009h used the notation 
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for which two different notations should be used. 
D(APE) to refer to the dissipation of AGPE into IE. 

If one accepts the idea that D(APE) and W r ,turbuient actually represent two large 
conversions between AGPE and IE, and between IE and GPE r respectively, then the 
question arises as to how these conversions are actually achieved in reality since there is 
no direct conversion between IE and GPE in the classical description of the energetics of 
the fully compressible Navier-St okes equations? With regard to W ryrn i X i ng , the theory of 
the hydrostatic adjustment, e.g.. lBannon I (|l995l ) can be tentatively invoked to speculate 
on some of the physical mechanisms and processes involved. Physically, localised heat- 
ing/cooling anomalies due to molecular diffusion must cause localised pressure anomalies, 
which will propagate as acoustic waves whose energy can be converted into kinetic energy 
via the compressible work of expansion/contraction and then ultimately into gravitational 
potential energy via the buoyancy flux pgw. If this is what indeed happens, then it is 
interesting to note that the Boussinesq/anelastic approximations implicitly assumes that 
of all the internal energy lost to mechanical energy, all of it goes toward increasing the 
gravitational potential energy, whereas in reality, it seems plausible that some of it could 
go toward increasing the turbulent kinetic energy, thereby acting as a positive feedback 
on turbulent mixing. The latter hypothesis warrants further research, as acoustic waves 
are kno wn to be c a pable o f generating mean flows via acoustic streaming for instance, 



Lighthill I |l978al) : IVanneste fc Biihler I (j201lh . What are the physical mecha- 



e.g., see 

nisms involved in the opposite conversion whereby AGPE is ultimately dissipated into 
IE by molecular diffusion is less clear, because physically AGPE can only be converted 
reversibly into KE by construction. As a result, the only way to dissipate AGPE into IE 
by molecular diffusive processes seems to require converting AGPE into KE reversibly, 
then KE into IE reversibly, presumably in the form of acoustic waves, and finally remov- 
ing the latter by thermal dissipation. The problem with this hypothesis, however, is that 
thermal dissipation is generally found to be a significa ntly less effective way to dissipate 
acoustic waves in liquids than the bulk viscosity, e.g., Lighthilll ( 1978bl) (thermal dissi- 
pation is more effective in gases, however). However, as the bulk viscosity dissipates the 
divergent velocity component, it is absent as a dissipation mechanism in the Boussinesq 
and anelastic approximations. If the physical mechanisms underlying the AGPE dissipa- 
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tion into IE are indeed related to the dissipation mechanisms of acoustic waves, then it 
is probably not possible to ascertain that only molecular diffusion is involved in reality. 

The overall conclusion is that the diabatic effects due to molecular diffusive processes 
in turbulent stratified fluids seem to give rise to nontrivial and potentially large conver- 
sions between mechanical energy and internal energy even in fluids traditionally regarded 
as incompressible or nearly incompressible, and that such conversions are actually sup- 
ported by such models as the Boussinesq and anelastic approximations even if this is 
still largely unappreciated. The main consequences is that thermodynamics, compress- 
ible effects, and the divergent component of the fluid velocity play a potentially more 
important role than traditionally assumed for understanding the physical processes and 
mechanisms ultimately involved in the energetics of turbulent mixing in stratified fluids. 
The present results, which suggest that the Boussinesq and anelastic approximations 
can support large conversions between mechanical energy, help rationalise why such ap- 
proximations appears to do so well in simulating turbulent stratified flows. On the other 
hand, the present results also suggest that real turbulent stratified fluids should exhibit 
differences with Boussinesq and anelastic fluids, as there must be a limit beyond which 
neglecting the effects of a divergent velocity may become noticeable for instance. Mak- 
ing progress toward clarifying these issues will probably require further detailed analysis 
of the energe tics of th e comp ressible Navier-Stokes equations along the lines recently 
developed by iTailleux I (|2009lh and direct numerical simulations of turbulent mixing in 
fully compressible stratified liquid flows resolving acoustic waves emitted by molecular 
diffusive heating/cooling and their dissipation mechanisms. Laboratory experiments, of 
the kind illustrated in Fig. Q] might also help. 
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